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Abstract 

The Stochastic Vacuum Model for QCD, proposed by Dosch and Simonov, is 
fused with a Worldhne casting of the underlying theory, i.e. QCD. biiportant, non- 
perturbative features of the model are studied. In particular, contributions associated 
with the spin-field interaction are calculated and both the validity of the loop equa- 
tions and of the Bianchi identity are explicitly demonstrated. As an application, a 
simulated meson-meson scattering problem is studied in the Regge kinematical regime. 
The process is modeled in terms of the "helicoidal" Wilson contour along the lines 
introduced by Janik and Peschanski in a related study based on a AdS/CFT-type ap- 
proach. Working strictly in the framework of the Stochastic Vacuum Model and in a 
semiclassical approximation scheme the Regge behavior for the scattering amplitude 
is demonstrated. Going beyond this approximation, the contribution resulting from 
boundary fluctuation of the Wilson loop contour is also estimated. 

PACS: 12.38.-t; 12.38.Lg, 12.38.Aw. 



1. Introduction. 

The confrontation of non-perturbative issues associated with dynamical processes consti- 
tutes a problem of great importance which merits definite attention if QCD is to attain the 
status of a complete and fully self-consistent theory. Clearly, the most concrete advancement 
in formulating a non-perturbative casting of QCD is traced to Wilson's proposal [1] which 
paved the way for the lattice formulation of gauge field theories in general. Remarkable 
results have been produced, especially in relation to the study of static properties of hadrons 
[2], finite temperature properties of the theory, etc. 

On the analytical front, important theoretical progress, relevant to non-perturbative, 
dynamical explorations of QCD, has been achieved within the context of the loop equations 
[3-6], while, in recent years, (super)string theory has, through the AdS/CFT conjecture [7- 
9], opened new pathways for approaching the non-perturbative domain of QCD, albeit in 
the sense of some supersymmetric version of the theory and within the context of unified 
schemes. 

Generally speaking, the nontrivial aspects of QCD as a relativistic gauge field theoretical 
system, stem from the fact that the non-abelian gauge symmetry entering its description 
incorporates an inherent non-linearity even before interaction terms with matter field agents 
are introduced. Theoretical schemes aiming at a heads-on analytical confrontation of non- 
linear quantum field systems do, of course, exist, possibly the most concrete one being 
expressed in terms of the infinite battery of the Schwinger- Dyson equations. Even in this case 
however, the relevant computational procedure for the solution of these (integral) equations 
is not only methodologically complex but, more importantly, unless totally summed, there is 
no a priori guarantee that they are in position to capture the full non-perturbative content 
of any given field theoretical system and/or describe its expected various phases. 

A notable field theoretical approach aiming at the study of non-perturbative issues in 
QCD, such as confinement, chiral symmetry breaking etc, has been proposed by Dosch and 
Simonov [10-12] and goes by the name of Stochastic Vacuum Model (SVM). By design, 
the construction of the model takes into consideration the non-trivial structure of the QCD 
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vacuum state [13], while, at the same time, it secures a role for the Stokes' theorem (non- 
abelian casting thereof) through which electric-magnetic duality issues can be fully taken 
into account. The basic building blocks of the SVM scheme are the, so called, fi^eld strength 
correlators, the definition of which will be given later, while for its computational strategy it 
employs the so-called Field (strength) Correlators Method, FCM for short. For the reader not 
familiar with the SVM, we hope that the information provided in this paper will sufficiently 
illustrate the reasoning behind its definition, as well as its properties and physical content. 
For a deeper insight to the model, one is referred to the original papers [10-12] and/or the 
excellent review articles [14,15] which also present a variety of its applications. 

The characteristic aspect of the present work is that it adopts, as basic methodological 
tool, the Worldline casting of gauge field theoretical systems [16] appropriately adjusted to 
the SVM. Our goal is to confront genuinely non-perturbative issues associated with dynam- 
ical processes of physical interest. At the same time -and at a purely theoretical level- we 
hope that the present effort will offer new insights further promoting the effectiveness of the 
SVM as a credible and viable theoretical tool for exploring the non-perturbative content of 
QCD. 

The aim of the first part of this paper (Sections 2 and 3) is to accomplish the task of 
carrying out cumbersome computations which reveal fundamental properties of the SVM 
that are of immediate relevance to our specific purposes and to establish the consistency of 
the model with the loop equations [3-6], as well as the Bianchi identity for QCD. Clearly, 
such an occurrence strengthens the credibility of the SVM as a theoretical construction which 
is consistent with QCD as a whole, i. e. in the sense that it goes beyond perturbation theory. 
Our main theoretical application will be realized in the second part (Sections 4 and 5)where 
we undertake the description of a, theoretically simulated, meson-meson scattering process 
at the high energy, small momentum transfer kinematical (Regge) regime. The relevant 
description of such a dynamical process necessarily "protrudes" into the non-perturbative 
domain of QCD as represented, in our case, by the SVM. 

Our presentation is organized as follows. In section (2), we discuss general aspects which 
prepare the "fusing" of the Worldline description of a non-abelian gauge field theoretical 
system, such as QCD, with the SVM. We start by focusing our attention on a situation 
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where a matter field entity of spin j interacts with a set of non-abehan gauge field modes as 
it propagates along a closed contour in (Euclidean) space-time, hence subjected to a spin-field 
interaction. The basic dynamical content of the process will be displayed by two alternative 
formulations. The first focuses on the (closed) Wilson contours traced by the particle entity. 
The second is based on a shifted field strength tensor that is integrated over an arbitrary 
surface bounded by the closed contour. Obviously the (non-abelian) Stokes' theorem plays 
a central role in relating the two descriptions, an occurrence of central importance to our 
purposes, given that the Stokes' theorem enters the SVM scheme in a major way. 

We shall subsequently introduce the so-called cluster expansion [14, 15, 17], which em- 
ploys the field strength correlators, the basic dynamical quantities of the overall description. 
In fact, the cluster expansion is the essence of the stochastic nature of the model and provides 
the key element for quantifying the stochastic vacuum hypothesis. Once the aforementioned 
task is accomplished, we shall be in position to derive, as a first general result, the equation 
that determines the surface over which the two-point correlator must be integrated. 

Section 2 deals also with fairly demanding calculations within the context of the Worldhne 
formalism whose first result is the derivation of an explicit expression for the spin factor. 
The latter represents the genuinely non-perturbative, spin-field interaction dynamics and 
necessarily enters [18] the analysis of the meson-meson scattering process. Section 3 is 
devoted to the verification of the loop equations and of the Bianchi identity in the framework 
of the adopted approach. We consider these results to be of importance since they solidify 
the credibihty of the SVM as a construction which reproduces sound, theoretical properties 
of QCD and demonstrate the compatibility with its non-perturbative content. 

In Section 4 we present a semiclassical calculation of a simulated meson-meson scattering 
amplitude. The Wilson loop that carries the dynamics of the process is a helicoidal embedded 
in a 4-dimensional background. The calculation is performed in the framework of the SVM 
and lead us to a Regge-type behavior for the amplitude valid in the physical region of the 
scattering. Corrections related to the fluctuations of the aforementioned Wilson contour 
will be discussed in Section 5. In the same Section the contribution of the spin- factor is 
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examined. Some concluding comments will be made in the closing Section. 
2. Worldline Formalism and Field Strength Correlators. 

In this section, certain basic features of the Field Correlator Method [15] will be reviewed 
in the context of the Worldhne casting of a Quantum, Gauge Field Theoretical System in 
interaction with a matter particle mode of a given spin j. More explicitly, the particle entity 
is taken to propagate along a given, closed, contour (worldline) while interacting with a 
dynamical set of non-abelian gauge fields A. According to our introductory discussion, the 
main objective is to introduce the necessary tools to facilitate non-perturbative, theoretical 
explorations of QCD in the framework of the SVM. 

One of the most important advantages of the Worldline formalism is that it allows one 
to reduce the physical amplitudes to weighted integrals of averaged Wilson loops [16]. In 
this connection, one can apply powerful techniques, such as the cluster expansion, both 
at the perturbative (in the sense of series resummation) and at the non-perturbative level. 
For the latter, the Worldline formahsm proves to be quite crucial because one can develop 
methods based on background gauge fixing strategy [18-22], which enables one to treat the 
non-perturbative fields as background. 

Let us, then, consider a particle entity of spin j propagating from some point and back 
to the same point while interacting dynamically with a non-abelian gauge field system A. 
The basic structure of the quantum mechanical amplitude associated with such a process is 
written in the Worldline formalism (all indices suppressed; Euclidean formalism adopted) as 
follows [16] 

K{L) = Tr J Vx{T)exp(^-^£dTx'^^(^Pexp(^iJ^dTX-A + Lj^dTJ-F^^ , 

a;(0)=a;(l) 

(1) 

where it should be noted that the parameter L has dimensions m^^ and must be integrated 
over through a weight factor exp(— Lm^) in order to obtain a result with physical content. 
The matrices J^i, stand for the Lorentz generators, pertaining to the spin of the propagating 
entity. Accordingly, the last term represents the spin-field interaction. 
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The above expression for K{L) can be recast into the form [16] 



K{L) ^Tr J Vx{t) exp (-^ ^ drx!^^ P exp drJ ■ — j (^P cxp (^i dr x ■ aJ 

(2) 



x{0)=x{l) 

where 



V 



= \im [ dhh pr p- (3) 

Sa^Axir)) .-o_y sx^ (r + |) 5a;, (r - |) 

defines a regularized expression for the area derivative [3-6]. 

Strictly speaking, expression (2) has a well defined meaning only for smooth [23] loops. 
On the other hand, when such expressions are used for the purpose of describing physically 
interesting processes the contour is forced to pass through points Xi where momentum is 
imparted by an external agent (field). Such a situation is mathematically realized by inserting 
a corresponding chain of delta functions S[x{Ti) — x] in the integral, which produces a loop 
with cusps, in which case the action of the area derivative operator entering (3) must be 
understood piecewise, i.e., 

Pexp . Aj = . . .Pexp (^^JirJ ■ AJ Pexp (^^JdrJ ■ . (4) 

The first step towards the application of the FCM is taken by employing the non-abclian 
Stokes' theorem [24] with the help of which one can write (the symbol Ps stands for surface 
ordering [14]) 



W[C] = -^Tr (^P exp (i £ dx • ^) = -^Tr ^P, exp i J dS^,{z)G^,{z, x^) 



A 

(5) 



The above expression is valid for any loop C with disc topology, irrespectively of the 
surface S{C). We also mention that for the area element we adopt the standard expression 

dS,, = \d^i^t,,{0, V(0 = -^e'^'daz.d.z,, a, 6 = 1, 2; {i\ e) = (r, s). (6) 

Finally, in relation (5) we have set [14,15] 

Gt,u{xo-, z) = (f){xo, z)Ff,^{z)(j){z, xo) (7) 
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with 

(p{z, Xo) = P exp I i j dw ■ A 

\ Xo ) 

a phase factor [14,15] which is a parallel transporter known also, in the SVM nomenclature 
[17], as connector. The reference point is chosen arbitrarily on the surface S\ arbitrary is 
also the curve that joins the points and z. 

It can be proved [24] that (5) depends neither on the surface nor on the contour used to 
define the connector (8), as long as the non-abehan Bianchi identities are satisfied. In the 
loop language such a requirement can be cast into the following statement: 

The relation 



-Tr ( P exp 



L sic) 



/ dS^^{z)G {Z,Xo) 



) =0 (9) 

A 



SzxiO 

is valid independently of the surface choice provided that 

which corresponds to the Bianchi identity for the gauge system [6]. 

The simplest way to prove the non-abelian Stokes theorem is to adopt the contour gauge 

[15] 

Ai^{x) ^ dsdsZ^{s,x)-^zx{s,x)F^x{z{s,x)), (11) 

with {zn{s,x), s e [0, 1]} an arbitrary, smooth curve from the reference point Xq to some 
point x: 

z^{0,x) ^ xofj,, z^{l,x)^x^. (12) 
Indeed, using Eq. (11) one can immediately see that 

^ dXf,A^,{x) = J dr J dsdsZ^drZyF^^{z) dr J dse"'^daZ^dbZ^F^^{z) = ^J (fi\f9ti^u{z)F^^{z). 

(13) 

For the gauge choice (11) the vector potential satisfies the condition tn{x)An{x) — 0, with 
tn{x) — dsz{s,x)\s=i, which implies that the connector (8) can be considered as the unit 
matrix in the contour gauge. In any case, the presence of the connectors in (7) guarantees 
gauge invariance. 
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The next step towards the apphcation of the FCM is to introduce the so-called cluster 
expansion [10-12,14,15,17] for the Wilson loop, formally written as 

1 / °° i"^ /■ ^ 

W[C] = — Tr exp IJ^—J dS^^^^ ■ ■ ■ dS^^^^ « G^„,.„(^„, xq) ■ ■ ■ G^^y^{zi,XQ) » , 

(14) 

where the symbol <<•••>> translates as follows: 

<< 0(1) »=< 0(1) > 

« 0(1)0(2) »=< P,(0(l)0(2)) >-\< (0(1) >< 0(2) >-\< (0(2) >< 0(1) > 

« 0(1)0(2)0(3) »=< Ps (0(1)0(2)0(3)) > -\{< P,(0(l)0(2)) >< 0(3) > +cycl.perm.) 

+|(< 0(1) >< 0(2) >< 0(3) > +cycl.perm.) 
H . 

(15) 

and is reminiscent of the cluster expansion in Statistical Mechanics. It is pointed out that, 
due to the color neutrality of the vacuum, expectation values of all correlators in (14) are 
proportional to the unit matrix in color space. This makes color ordering unnecessary. 

The formula of Eq.(14) quantifies the formulation of the SVM. It turns out [14,15]that the 
most important contribution to the cluster expansion comes from the two-point correlator: 

aS,p(^ - z') = j^Tr{G,,{z,xo)Gx,{z',xo))A = -^Tr{F,,{z)<P{z, z')Fxp{z')<P{z' , z))a. 

(16) 

The above defines the field strength correlator, a quantity which constitutes the basic building 
block of the model. 

Some natural assumptions are incorporated in the above definition. The first has to do 
with the Lorentz invariance of the vacuum, which is explicitly indicated in the left hand side 
of (16) by the fact that the correlator depends on the distance between the points z and 
z'. The second is that, on the other hand, the correlator does not depend on (the gauge 
parameter) xq. This is a credible assumption, taking into account the fact that we have been 
working, from the very beginning, with a gauge invariant amplitude. Accordingly, the basic 
assumption of the SVM leads to the statement that 



W[C] oc exp 



I dS,,{z) J dSxf>{z')A^ 



S{C) S(C) 



(17) 
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where the surface element enters through the use of the (non-abehan) Stokes' Theorem. 

Two important points should now be made. First, the last expression is supposed to be 
valid in a certain limit. Explicitly, it is assumed [10-12,14,15] that the vacuum fluctuations 
establish a correlation length Tg beyond which correlations decay very fast. If A is an 
order of magnitude estimation for the two point correlator, relation (17) is considered as an 
asymptotic approximation which is valid in the limit TgVA — > 0. The second is that, while 
expression (14) does not depend on the particular surface one uses for the application of the 
non-abelian Stokes' theorem, approximation (17) does. Thus, the stochasticity assumption 
transforms relation (9) to an equation which determines the dominant surface in the cluster 
expansion. To quantify this statement we write 

= l J dS,Az) J dS,,{z')A^l^{z-zr (18) 

S{C) S{C) 

from which it is easily determined that 

= / dS,,{z') [a,Af ),^(i - z') + a.A(2,^(5 - z') + a.Ag,,^(5 - z') 

^"^^^^ S{C) 

(19) 

where we have written = The calculation of the derivative of the correlator (16) 
needs to take into account [15] the presence of the connectors 0(^, Xq). In Appendix A we 
show that 

d 

' 'z, xo) = iA^,{z)(j){z, xo) - il^{z, Xq), (20) 



dzu. 



with 



I^{z, Xo) = ! dt (dtUJn) ^H^, a;)F«A(a;)0(a;, xq) 



ijj — u!{t, z), 00(1, z) — z, a;(0, z) — Xq. 
Given the above, we conclude that 



(21) 



d^Gapiz, Xq) = 0(xo, z)Dij_Fap{z)(f){z, Xq) + i[Iij.{z, Xq), Gap{z, Xq)] (22) 

and consequently write 

di^Gi^^{z, Xq) = (f){xo, z)D^F^^{z)(t){z, Xq) + i[I^{z, Xo), Gfj,u{z, Xo)], (23) 
where we have set G^i/ = ^e^^^^f^Gap- 



The above analysis establishes the following relation [15] among the derivatives of the 
correlator: 

1^.«..5^^(2)^^ ^ d^Afl,^ = ^Tr{Dj^^{z)F,,{z')cf>{z\ z))^ - A.,,(^, z% (24) 
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where 

A«Ap(^, z') = ^Tr{F,,{z)I„{z, z')F^p{z')4>{z, z') - F,,{z)4>{z, z')F^,{z')I,{z' , z))j,. (25) 
Thus, if one assumes the validity of the Bianchi identities D^Ffj^^ — 0, one concludes that 

dA'^'l, + d.A^l,, + d^A^^l,^ = e'^-A.,,. (26) 
Accordingly, Eq (9) can be represented as follows 

= ^V^V / dS,,{z)e-^^A,,,{z -z) = 0. (27) 

S{C) 

It is worth noting that A^^xp is a three point correlation function and it would be identically 
zero if we were considering an Abelian gauge theory so that the above relation becomes, 
really, an identity, telling nothing about the particular surface involved in Stokes' theorem, 
an expected result given that relation (17) is exact in the framework of QED. An extensive 
discussion of the physical content of the correlator A^xp can be found in [14,15]. According 
to the analysis presented in the aforementioned references confinement in QCD occurs due to 
the non-zero value of the (non-abelian) correlator A^xp- Eq.(27) indicates that this correlator 
also defines the relevant surface on which the two-point correlator "fives" . 

Now, it has been demonstrated [10-12,14,15] that, in the asymptotic fimit | S — z |» 
Tg, Eq.(27) determines the surface 5" bounded by the contour C as the minimal one. To 
demonstrate this we employ the following general Lorentz structure representation [14,15] 
for the two-point correlation function: 

a!2,Ap(^) = iS,xSup - SMD{z) + |1 A[(^,5,, - Zp6.x)D^{z)] - (/. - i^)| , (28) 

where we have set z = z — z'. 
It is easy to see that 

1^,...5^^(2)^^ = e^'^'^d^D = A.,,. (29) 
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With the help of the above relation Eq.(27) can be cast into the form 

r 

j—r^ I dS,,{z) I dS^,{z')D{z - z') = 0. (30) 

'^'^y^i s{c) s{c) 

The functions D and Di have been measured in lattice calculations [25] and have been 
found to be of very fast decrease as \ z — z' p— > oo. In the considered region both of them 
were found to be of the form / ^^^1^^ and that they go exponentially fast to zero for 
\ z — z' \> Tg. In this region, we write: 

Z{C') ^ ^(0 + (r - O'^daZiO, \Z-Z' - 0"(r - O'gab. (31) 

In the considered limit and taking into account that t^^t^y = 2 we find 

aS + C(T„^A)1 = 0, (32) 



where the string tension 

(^ = yI d^'^Diw''), (33) 



2 

has been introduced and where we have also written 



S = j dSv^) (34) 

for the area of the surface bounded by the Wilson curve. In the last equations we used 
the dimensionless parameter — ■^-{z — z')^ and we have written Qabiw) — daW/^dbUUfj, the 
induced metric. It accordingly follows that in the hmit — > 0, the surface on which the 
two-point correlation dominates is the minimal one. 

Now wc turn our attention on the spin factor, whose role is to incorporate the spin-field 
interaction in the framework of the worldline formalism. We mention that in the present 
work we shall be dealing with massive fermions. It should be noted, at the same time, 
that the spin factor can also be extended [21,22] for the case of, massless, bosons of spin-1. 
Accordingly, the calculations to be presented in this section can easily be extended to bosonic 
fields. 

We start by inserting into the SVM formula for the Wilson loop, cf. Eq.(17), into the 
worldline integral expression (1), the basic goal being that of calculating the spin-field inter- 
action with the help of the area derivative operator defined in (3). Once this is accomplished 
the stage will be in place for performing specific calculations of physical interest. 
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We start by introducing the spin factor by formally casting Eq.(l) into the form 

K{L) =Tr J Vx{t) exp (^-^ drx^^ *(j)[C]e-^[^], (35) 

x{0)=x{l) 

where the quantity 

$(^) [C] = e^[^]p exp I'dTJ-j-] e-^[^] (36) 
\2 Jo da J 

defines the spin factor characterizing a particle entity propagating on the Wilson curve A[C]. 
Its calculation is not trivial and requires a number of steps the first of which is to determine 
the action of the area derivative operator on A. 

The first objective of the computation of the spin factor is to study the change on A[C] 
induced by an infinitesimal variation of the boundary. The relevant problem is formulated 
as follows 

SA 



= Jds[a(n,s)z^{n,s)]' I dS,s(z')A^ZyS 

5(C) 



+ 



1 



+ Jdsa{n,s)z^{n,s)z'^{n,s) J dS,s{z') {d.A^^s + daAf^,,s) , (37) 

5(C) 

where the dot denotes (partial) derivation with respect to r, while the prime derivation with 
respect to s. Moreover, the correlators depend on the distance | z{ti,s) — z{t',s') \ and we 
have written = Sa^cS(r - ti)q;(ti, s). 

Using Eq. (26) equation (37) is recast into the form 

1 

J ds[a{Ti,s)za{Ti,s)]' j dS^s{z')A^^l ,^ 



5A 



^' 5(C) 



1 



- Jdsa{n,s)z^{n,s)z'^{n,s) J dS^siz'^d^A^^^l^s - ^'""^^^-yus) = 

5(C) 

= Un) I dS,s{z')A'^^l,sHT,)-z{T',s')] + 

5(C) 

1 

+ Jdsa{n,s)za{Ti,s)z'^{Ti,s) J dS^s{z')e^''''^ A,^s[z{n, s) - z{t' , s')]. (38) 

5(C) 

The last term in the above relation is zero on account of condition (27): 

Za{ri, s)4(ri, s) J dS,s{z')e''''''^A,,s[z{n, s) - z{t', s')] = 

5(C) 
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= \^)tap{z) j dS,s{z')e^''''PK,5{z-z'). (39) 

5(C) 

We have consequently determined that 

5A 



J^^=ia{T) / dS,s{z')^':i^,[x{T)-z{T\s')] (40) 



S{C) 

In order to find the area derivative, c/. Eq.(3), we need to calculate the second functional 
derivative of ^4, at the points x{ti) — x {t -\- ^ and x{t2) — {j ~ Ftom the definition of 
the area derivative we also surmise that only terms ~ 5'{h) are relevant. Accordingly, it is 
straightforward to surmise that 

5A 



5a,Mr)) ^^^^ 
It is also easy to determine that 

s s 



J dS,s{z')A^Xs[<^)-<^'^'')] (41) 



while all higher derivatives give null contribution. 
On the basis of the above analysis we determine 

1 



^ = A(l,,^,Jx(r2)-x(rO], (42) 



iL 

5(C) 



$0)[C] = 1-'— Jdn J dS- A^^\z -xi)-J + 



■ T 2 ^ "^2 

+ J '^^2 / dri -J ■ A^^\x2 -xi)-J+ J dS- A^''\z -X2)-J j dS- A^^\z - xi) ■ J 

5(C) 5(C) 





3 1 T3 T2 



+ 



+ (y)' J drs J dT2 j dn [j ■ A^^\x2 -xi)-J j dS- A(2)(^ - Xg) • J + 

5(C) 

+ J ■ A^'^\x^ - X2) ■ J j dS ■ A^'^\z - xi) ■ J - j dS ■ A^^\z - xs) ■ J X 

5(C) 5(C) 

X J dS ■ A^^\z - X2) ■ J J dS-A^^\z-xi)-J+---.] (43) 



In the above expression we have omitted terms in which the field strength correlator A^^^ 
depends on distances between two non-successive points. For example, the last line of the 
above equation does not include the correlator A'^'^\x{ts) — x{ti)). The reason is, that such 
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a correlator [14,15,17,25] is assumed to behave as / ^ / f a:^ '-^"^a^-' ), which, in 



turn, means that its contribution is suppressed by powers of T^vA. Accordingly we obtain 
the following expression for the spin factor: 



$(^)[C] = Pexp 



^L}■ r .„ .... . ' 



j dr j dS ■ AP\z-x)- J + — j dT2 j dTiJ ■ AP\x2-Xi)- J 



S{C) 



14) 

We have arrived at a result of considerable interest for our purposes, which we would like 
to analyze further. Let us start from the second term in the exponent which is of special 
interest [26] and represents the interaction of the quark color-magnetic moment with the non- 
abelian background. In particular, let us refer to the representation (28) of the two-point 
correlator, which we rewrite in the form 

^SapC^) = {^u.\^vp - Sfj,p5^x)iD + Di) + {x^xx5vp - x^Xp5^x - (a« ^ t^))D[. (45) 

where we have denoted D[ — -^Di[x). 
Using expression (45) we find 

^m^a2],,,(x) J,, = 2J\D + D,) + 4( J • xfD[. (46) 

For the case in hand we have that J^^ — |[7//, 7;/], so we can easily determine that 

>/M.A(fj,,,(x) = Q{D + D,) + ^x^D[. (47) 

Expression (47) is positive definite and is associated [27] with the ghost tachyonic pole which 
appears in the fermionic, or the gluonic, propagator -an issue we shall not discuss further in 
this paper. 

With the help, now, of result (41) the first term in the exponential (44) can be recast 
into the form 

jir j iS.^^^\.-.).J = jirJ^^^ (48) 

5(C) '"^^ ^ 

It is convenient, for the applications we have in mind, to rewrite the above relation, by 
referring to the variation of A[C\ as it has been computed in Eq.(40): 

5A[C] 



Sxi^{t) 

14 



^9Mr)]- (49) 



The above function is reparametrization invariant, thus x^g/j, — 0. 
Taking into account the result displayed in Eq.(40) we can write 



An obvious solution of the above equation is 



>m 1 (51) 



5ai^,y{x) x'^ 

It can also be shown [28] that the above is the only possible solution. The proof follows, 
basically, dimensional arguments: Making the change x ^ \x and r — > Ar, the area deriva- 
tive of A[C] scales like The same scaling behavior goes for the function as can be 
easily concluded from Eq.(40). Thus, the term appearing in (51) has the right scaling prop- 
erties. Any other term must be an antisymmetric combination K^j^i, which scales as ^ and 
is perpendicular to the velocity, i.e., x^K^j^^ = 0. However, such a combination expressible 
in terms of the boundary x{t) cannot be found. It thereby follows that the first term in the 

exponential which defines the spin factor reads 

1 1 . . 

jdr j dS-A^^\z-x)-J = JdT^^^^^^^^^^^J^,. (52) 

S{C) ^ 

3. Loop Equations, Bianchi Identity. 

In this section we shall proceed to assess the capacity of the SVM to expedite non- 
perturbative investigations in QCD by examining whether Eq.(17), as formulated within the 
framework of the stochastic approach, satisfies the Polyakov/ Makeenko-Migdal equations 
formulated in loop space. The latter, constitute the most credible proposal for achieving 
a non-perturbative casting of the theory, equivalently, one which provides a solid basis for 
conducting non-perturbative investigations within its framework^. In addition, we shall 
explicitly demonstrate the validity of the Bianchi identity. We claim that these properties 
must be satisfied, if one is to assert that the expression for the Wilson loop, as given in Eq. 
(17), is to constitute a credible approximation to the full theory. 



^One might consider this approach as the continuum space casting of lattice gauge theories, which defines 
the 'going' standard for such investigations. 
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Now, the loop equation for a contour without self-intersections can be stated [6] as follows 



r+e 



r^^[C] = lim / df——-- ——W[C] = 0. (53) 

Its verification constitutes the first, as well as simplest, test the SVM must pass. To this end, 
let us insert Eq.(28) into Eq.(41) whereupon, using the fact that the boundary is a closed 
contour, we determine 

S 



+ 2/0 "^^'i^i^i^'^M^') - ^-(^)) - (/^ ^ ^)]DMr') - x{t)]. (54) 

Our next step is to take the functional derivative of the above equation. Our task becomes 
relatively easy as we notice, from Eq.(53), that we only need those terms which contain the 
delta function S{f — t). Accordingly, we obtain 

+ 2 / dAx,{T'){x^{T') - x.(r)) - (a. - ^)\j^^D,{x{t') - x{t)) + 

+ -5{t -f)j dr'x,{r')Di{x' - x). (55) 



We now write 



^ 7^/ / NX r/ ^,dD(z — x) ^. ^.dD(z — x) 

D{z{0 - x{t)) = d{T - r)— ^- ^ = -d{T - t)^- ^ (56) 



5xfj,{f) dx^ dZfj, 

and similarly for the derivative of Di. Thus 



'^''sJdry^ ^ -/.Ve-.„.,(oa...(r)^i.(.-.) - 

^y'dT'[i:^(T0(M^')-^.(^))-(/^^^)]^A(2:'-x) + ^y'dx,L'i(x). (57) 



— iJ- 

Since the boundary is a (closed) loop, we conclude that the first term takes the form 

2m 



J d'C'e'^'daZ.iOd.z^iOg^Diz -x) ^ j d'C'e'^'daDiz - x)d,z,{e) = 

1 

^ - j dT'x^{T')D{x' -x) ^ - ^^dx^D{x).{58) 
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With the same reasoning we have for the second term 

I dr' [x^{t') {x,{t') - x,{t)) - (// ^ u)] g^D.ix' - x) ^ 

1^ 



= - 1 dT'x,{T'){x^{T') - {x^{t))—D,{x' -x) = -2j^dx,\x ^ -^D,{x).{bQ) 
Accordingly, 

(60) 



^ 5(Tfj,v{x{T)) JC 



D{x) + Id,{x)+ I X |2 



2 ^ ' ' 9 I X |2 



Given that the functions D and Di depend only on the distance | a; |, the right hand side of 
the above equation vanishes and the loop equation is satisfied. In fact what is known is that 
the minimal area satisfies the Makeenko-Migdal equation aymptotically. The above result, 
on the other hand, can be considered as new in the sense that it is an exact result in the 
framework of the Stochastic Vacuum hypothesis. 

Our next concern is to confirm the Bianchi Identity (BI) in the framework of the SVM. 
The reason we are interested in such a confirmation is based on the fact that the zigzag, 
or backtracking, symmetry characterizes the Wilson loop functional: It is invariant under 
reparametrizations of the form x{t) x(q;(t)), even if a' < 0. It can be shown [24] that 
a Stokes type functional, of which the Wilson loop is a prime example, which respects the 
aforementioned symmetry also satisfies Eq. (10), which can be written in the form 

The above relation forms the bridge between the zigzag symmetry and the BI in the frame- 
work of a gauge field theory. 

As it is, now obvious the area S = J d^w^ does not respect the zigzag symmetry and 
in this sense e"'^'^ is not a good representative for H^[C]. On the other hand, the action 
A[C] in the stochastic approximation is invariant under zigzag parametrizations, but due 
to the truncation it is not obvious that it is a Stokes type functional. We are thereby 
obliged to confirm that the BI is exphcitly satisfied and that, consequently, the stochastic 
approximation produces Stokes type functionals. To this end we first observe, using Eq. 
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(41), that 



5a^y{x{T)) 



5 




5(C) 



I dS^f,{z)d,A^;i^{z-x)\e 



AlC] 



(62) 



Using Eq.(29) we find that 



(63) 



5(C) 



where we have set 




dSapA^ap{z - X). 



(64) 



5(C) 



Now, we have estabhshed that the surface over which we integrate the correlators is deter- 
mined by Eq. (27). Multiplying that equation with e'^'^'^'"' we determine 



The above equations form a homogeneous system whose only solution is X^, = 0, a result 
which, together with Eqs. (61) and (62) confirms the validity of the BI within the framework 
of the SVM. 

4. Meson-Meson Scattering. 

In addition to confinement, which constitutes a profoundly non-perturbative problem, 
there do exist specific dynamical processes, whose theoretical confrontation also calls for non- 
perturbative methods of analysis. One such situation arises in connection with the theoretical 
description of high energy scattering amplitudes for which the soft sector of the theory is 
involved. Prom the experimental point of view, one such case arises in connection with 
Regge kinematics, entering directly the theoretical description of, among others, diffractive 
and low-x physics processes. In this section wc shall study a simulated case of a meson-meson 
scattering process whose quark-based description is of the general form: 



Adopting a standard picture, already employed in the QCD literature (see, for example 
[29-31]) according to which quark 1 from the first meson and antiquark 2 from the second 




(65) 



(11) + (22) ^ (33) + (44) 
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meson are very heavy, in comparison to the incoming total energy, hence their worldhnes 
are considered to remain intact from the gluon field action and can be described in the 
framework of the eikonal approximation. The light pairs 1,2 and 3,4, on the other hand, 
are annihilated and produced in the ^-channel, where the eikonal approximation is not valid 
and a full treatment is called for their description. In the Worldhne framework the process 
is schematically pictured in space-time by the straight eikonal lines (1 — > 3) and (2 — > 
4), describing an intact quark and anti-quark and by the curves (1 2) and (3 4) 
which correspond, respectively, to the annihilated and produced quark antiquark pairs. The 
structure of the field theoretical amplitude can be written as follows (see Fig.): 

G{x4,X3,X2,Xi) = {iSF{x4,Xs \ A)iSF{x3,Xi I A)iSF{Xi,X2 I A)iSF{x2,X4 I A)) A- (66) 

In the above expression {Sf is the full fermionic propagator which, in the framework of the 
Worldline formalism, assumes the form [16] 



POO r 2 r ~i / '^'^^^ 

iSF{y,x\A)^J^ dLe-^"^ J Dx{T)e « 



x(0)=x 
x{L)=y 



m 



$(^/2)(L,0)Pexp ii j drx ■ a \ , 



(67) 

where ^^^^ is the so-called spin factor for the matter particles entering the system. For us, 
it means that j = \- 

Inserting the above formula into Eq.(66) we find 

G(a;4, xs, X2, ^i) = H / ^^^^(^^ ' Ti-i)e-(^^-^^-i)"^' /' /' Dx{r)5[x{r3) - X3]5[x{t2) - x^] x 

1=1 n 



x(0)=X4 



x5[x{ti) — xq] exp 



T4 



dr x^{r) 



(spin structure) (^P exp (^i dx ■ A 



(68) 



where the term spin structure corresponds to the following expression 



(spin structure) = JJ 



1=4 



(69) 



In principle, the Wilson loop appearing in Eq.(68) incorporates the dynamics (perturba- 
tive, as well as non-perturbative) of the process. In the framework of the SVM it assumes 
the form 
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e-^[^l (70) 



'Pexp(^t£dx-A^'^ =cxp -i J dS^,{z) J dSxp{z')A^^lxpiz - z') 

^ ^ [ s{C) S{C) 

In the present Section we are going to calculate the amplitude (68), using the above 
expression which gives the structure of the Wilson loop in the framework of the SVM. The 
particular method to be adopted is a kind of a "semiclassical" approximation based on a 
combined minimization of the action A[C]- see Eq.(39)- with respect to the surface S[C] and 
of the surface S[C] with respect to the boundary C. The reasoning behind this procedure 
is that, according to Eq.(68), in order to obtain the full amplitude it does not suffice to 
determine the minimal surface bounded by a given specific contour but one needs to proceed 
even further and sum over all possible boundaries with a weight of the form 

S[x]^-jdTx^ + A[C]. (71) 



The above described approximation will allow us to determine the dominant contribution 
to the Worldline integral (68) in the stochastic limit Tg\fK <S 1. 

The variation, g^[x{T)], of A[C] under changes of the boundary is given by Eq. (39). 
Accordingly, the correlator contributions become stationary for the "classical" trajectory 

= 0. (72) 

Using the expansion for the correlator according to Eq.(28) it is easy to see that 

9h[x{t)] = 2Xa{T)Rai^[x{T)] - ^Xa{T)Qa^,[x{T)], (73) 

with 

RaM^)]= j dS^^{z')D[x{T)-z{s'y)] (74) 

5(C) 

and 

QaM-^)] = j dr' [i:^(r') (x,(r') - x^{t)) - {fi ^ a)] Di[x{t) - x{t')]. (75) 

It is worth noting that the above expressions are reparametrization invariant. Also, in 
the last relation the integration covers the whole range of the r variable. Following Refs. 
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[30,31] the minimal surface bounded by two infinite rods at a relative angle 6, has (in four- 
dimensional Euclidean space) the shape of a (three-dimensional) helicoid, which is the only 
surface that can be spanned by straight lines [32]. In the considered process the eikonal 
lines 1 ^ 3, 2 ^ 4, play the role of the "rods", while the angle 6 is connected, via analytic 
continuation [33], to the logarithm of their total energy s. 

Given the above specifications, consider the following, helpful, parametrization of the 
boundary C: For < r < ti we have a straight line segment, x^^\ going from the point 
a;4 to the point X2- Introducing, moreover, for convenience the length 2T =\ x^^ — X2 \ and 
reparametrizing according to r ^ — T, we write 



xW = (r, 0,0,0), -T<T<T, (76) 

with x^^\—T) = X4, x^^\T) = X2. 

The second eikonal line x^^\t), T2 < t < T3, goes from the point Xi to the point x^ at 
a relative angle 9 with respect to x^^\ while a distance b (impact parameter) separates the 
two linear contours in a transverse direction. Introducing the distance 2Ti —\xs — xi \ and 
reparametrizing according to 

r^rJ^r-'J^) (77) 



we write 



T3 - T2 T3 - T2 



xJf)(T) = (-Tcos^, -rsin^, b, 0), -Ti < r < Ti, 



with (-Ti) = xi, (Ti) = xs. 

In the following we shall assume, just for convenience, that 

2T =1 X4 -X2 h I X3 - xi 1= 2Ti. 

For Ti < r < T2, we have a helical curve x^^\t), which joins the points X2 — x^^\ti) and 
Xi — x^^\t2), representing the exchanged light quarks. Performing, now, the change 
(7 = — - — (r — Ti), we write 

xjf) (cr) = ^0((7) COS y , 0((7) sin y , CT, j , < (7 < 6. (78) 
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The continuity of the boundary requires 

^^Ht) = x(f)(0) = X, and = = x^, 



or 



0(0) = 0(6) = T. 



(79) 



The final hehcal curve is a;*^^)(r), which, for ra < r < r4, joins the points ^3 = x'^'^^r-^) and 
X4 — x^^^ (T4) . Making one more, final, reparametrization, namely a 



T4-T3 



—(p[<J) cos — , — 0((7j sm — , cr, 1 , < cr < 0. 



I!''" ~ ''3) we write 



(80) 



Once again, Eq.(79) takes care of the continuity of the boundary. Now, the minimal surface 
is bounded by the (four) curves specified by Eqs.(76)-(80) and can be spanned by straight 
lines parametrized as follows 

T-T . T + T 



2T 



(4), 



(7 + 



2T 



^(f){a) cos y , ^(f>{a) sin y , a, j . (81) 



It can be easily proved that the surface defined by the above equation is minimal, irrespec- 
tively of the function 0: 



V9 



(i • z')z^ 



V9 



= 0. 



(82) 



One observes that the minimization of the surface is not enough for the complete spec- 
ification of the parametrization of the helicoid. Accordingly, we go back to Eq.(72), which 
determines the boundary that dominates the path integration (68). A first observation is 
that, due to the antisymmetric nature of Rafj, and Qa/x, the function vanishes when Xju(r) 
represents a straight line. Thus Eq.(72) is trivially satisfied for the eikonal sector of the 
boundary. Non-trivial contributions arc coming only from the helices x^^^ and x^^^ . One can 
simplify Eq.(73) by computing the leading behavior of the functions Raf^ and Qaf^ using the 
fact that the functions D and Di, as defined in the SVM scheme -and measured in lattice 
calculations [25]- decay exponentially fast for distances which are large in comparison with 
the correlation length Tg. In this connection and upon writing 

x{a') = x{a) + {a' - a)x{a) + ^{a' - afx{a) H , 
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we find, for the second term in Eq.(73), 



1 r 
2 

1 



{x )Xfj, — {x ■ x)x^ J da {a — a) Di 





X 



9 



+ 



X 



X 



X ■ X 



x^ 



Tgai 



+ 



(83) 



where^ 



1 

ai ^ J 



In the last expression we have used, as in Eqs.(33) and (34), the dimensionless parameter 
w =\ z \ /Tg. It must be noted that the coefficient l/ai is a small number. Taking into 

account that Tg is of the order of O.lfm and the string tension (see Eq.(33)) a ~ 0.18 GeV^ 
it is readily seen that l/ai = 0{aTg). 



Noting that 



z,{a,r^ T) = xf{a), z,{a,T ^ -T) = xf{a) 

1 

'XT' 



drZ^{a,T) = i^((7,T) = — [xjf)((7) - xf {a)], 



(84) 



the leading behavior of the first term of the rhs of (73) can be easily determined. One finds 



[X ■ Z) 



-X Z. 



^1 I / . \x ■ z 
2T\x \ [z„- ^—r-x 



dr' 



T 
2 



'h' 



D 



x' 



9 



+ 



X 



2 ^ j T 



+ 



(85) 



where we have introduced the parameter 

oo 

//^ = T"^ j dwD{w^) ~ 0(a). 



Thus, the function g takes, to leading order, the form 

1 



(86) 



9ix 



iT/i'x' { z^ - ^x„] - ^ 



X 



2 -^M 



2ai 



X ■ X 



^2 



(87) 



Now, we recall from its definition that the gf-function provides a measure of the change 
of A[C] when the Wilson contour is altered as a result of some interaction which reshapes its 



^We have omitted terms suppressed by powers of Tg 
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geometrical profile. In this sense, it contains important information concerning the dynamics 
of the problem under study. The structure of the s^- function, as it appears in the above 
equation, is quite general and exhibits its dependence, not only on the boundary but on the 
minimal surface as well. It is worth noting that this fact is strictly associated with the non- 
Abelian nature of the theory since the function D -and consequently /z^- disappears [14,15] 
in an Abelian gauge theory. 

Taking into account that for the helicoids parametrization the velocity x has three non- 
zero components, while x and i have only two, we conclude that Eq. (72) can be satisfied 
only if 

4T^Ji^x'z^ - ^x^ = 0. (88) 

Inserting in Eq.(88) the helical parametrization one easily finds that the function must be 
a constant. Taking, now, into account Eq.(79) we determine this constant to be the length 
T. It is then very easy to see that this result leads to the conclusions 

i;-i = 0, x-x^Q (89) 

and 

This equation cannot be satisfied in Euclidean space. In Minkowski space the angle 6 
becomes imaginary 9 — > —ix — —i In (^) (s is the total energy of the "heavy" quarks that 
form the "rods" and m(~ mi ~ ms) their mass and Eq.(90) has a positive definite solution: 

In the last equation we have analytically continued only the angle between the "rods" and 
not the parameter T. Our reasoning is based on the one presented in [31]: Eqs.(79) and 
(90) define a contour of steepest descent for the path integration. For such a contour the 
parameter T is determined to be imaginary (the formulation is still Euclidean and 9 is real) . 
After performing the T integration the only parameter remaining for analytic continuation is 
the angle between the "rods". Eq.(91) indicates that, equivalently, one can first analytically 
continue the angle variable to imaginary values leaving the T parameter real. In any case 
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it is obvious that the impact parameter must grow with the incoming energy: Tgb ~ Ins, a 
conclusion which is in agreement with the landmark result of Cheng and Wu [34]. 

The preceding analysis obviously repeats itself for the two helical curves x*-^-* and x^^^ 
and has led us to a specific parametrization for the Wilson loop, which plays the dominant 
role in the path integration in Eq.(68). We are now in position to determine the leading 
contribution to the action (71): 

1 '^^ 

S,,^-ldTxUT)+A[CU. (92) 



Our first step is to expand the second term of the integrand in powers of TgV A. The first 
term of such an expansion is the familiar Nambu-Goto string. The next term, which reveals 
the rich structure of the SVM, is the so-called "rigidity term" , representing the extrinsic 
curvature of a surface embedded in a four-dimensional [35] background: 

A[C] =a [ d'i^ +— f d^i^g'^'dat^.d.t^, + • • •, (93) 

where a is the string tension as defined in Eq.(32). 
The coefficient of the rigidity term reads 

Terms proportional to Tg entering the expansion in Eq.(93) will be considered negligible 
in our analysis. We have also omitted the term / (P^^/gR, since in two dimensions the 
curvature is a total derivative. Using the helicoids parametrization (81), with (p — T, the 
Nambu-Goto term in Eq.(93) takes the form 



/ d'^,/g= j drj dsJl + ^ = bT 



-T 



l-hp^-l — ln(^l-|-p^-|-p 



(95) 



where p — 

To proceed further we analytically continue to Minkowski space where we can use Eq.(91) 
to determine 



,T^. ^ - . 6 (l - ^ . + 0(T|) (96) 
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and 

P 
Thus 



— In [sjl+V^ + 



hT 



-iTx/h 



In 



2x y/8a^ 



+ 0{T'). (97) 



(98) 



In the same framework, the contribution of the rigidity term takes the form 



-T 



e 



f- 1 fe^ le^ 

"T / ds I = — H — — ' 



-In \^\/l+p'^+pj + -P\jl+P^ 



It follows that in Minkowski space we have 

1 



37r 



jj,!/ 



4ar 



-X- 



(99) 



(100) 



For the full estimation of the classical action, cf. Eq.(92), one should also take into account 
the presence of the classical kinetic term. Non trivial contributions come from the helical 
curves x^^\l 2) and x(^)(3 ^ 4): 



b^P 



4(r2 - Ti) 



+ 



4(r2 - Ti) 4(t4 - T3) 



(101) 



' ' 

Now we have to take into account that both T2 — ti and T4 — T3 must be integrated 
with weights e~('^2-n)mg g^^^ g-(T4-T3)m§^ respectively. These integrals, as it turns out, are 
dominated by the values T2 — Ti = T4 — = leading to a final kinetic contribution of 
the form 



2mo& I X 



mo 



(102) 



Here, mo(— m2 — ^4) is the (current) mass of the light quarks, thus the result expressed by 
(102) can be considered negligible. 

From the above analysis we conclude that 



Sr. 



3n 
- 1 — X 



(103) 



2x 4q;o' 

Putting aside, for now, the possible corrections to A[C] which arise from fluctuations of the 
boundary as well as the spin factor contribution, let us consider the result (103) as a whole. 



26 



except for terms ~ mass. To obtain the final expression for the scattering amphtude one 
must integrate over the impact parameter: 

I d%e^p (iq- b- ^b'^ cx exp {-^q'x) ■ (104) 

Combining (103) and (104) we find, for the scattering amphtude, a Regge behavior of 
the form s"k(o)*+"fl(o) with 

1 Stt 
«r(0) = TT- and Q;ij(0) = — . (105) 

Zna 4q;o 

The form of the slope a'ji{0), indicated in the last equation, is the same as the one 
obtained in Ref.[32]. In that work the authors have applied, in the framework of the SVM, 
a different method based on the path integral- Hamiltonian duality and consequently their 
result is confined in the region t > 0. In this sense our result extends the same value 
a^(0) = l/27rcT ~ 0.9 GeV^^ for all the values of square momentum transfer. 

It is well known [26,32] that the intercept ai^(O) receives significant contribution from non- 
perturbative corrections to quark self-energy. We shall comment on this interesting issue in 
the next Section. The result indicated in Eq.(105) does not take into account the aforemen- 
tioned corrections thus it is very sensitive to different lattice data or parametrizations. For 
example, using Ref.[25], the coefficient I/oq of the rigidity term is negative and one needs 
[32] the large (and negative) quark self-energy corrections to restore the phenomenological 
value of the intercept. On the other hand adopting a certain [17,36] parametrization for the 
functions D and Di (see Appendix B) one obtains for the string tension the value a — 0.175 
GeV^ and for the coefficient of the rigidity term the value 1/ao = 0.276. With these numbers 
we obtain for the Reggeon slope the value aj^(O) = 0.91GeV~^ and for the Reggeon intercept 
the value aij(O) = 0.65 in good agreement with the phenomenological values Q;^(0) = 0.93 
GeV-2 and aR(0) = 0.55 [37]. 

5. Boundary Fluctuations and the Role of the Spin Factor. 

As repeatedly mentioned in our narration, corrections to the amplitude (68), beyond 
semiclassical ones, are expected to arise from fluctuations of the boundary of the surface on 
which the two-point correlator "lives" . Fluctuations of the surface itself can be taken into 
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account by higher order correlators. This, in fact, is the big difference which distinguishes 
the SVM approach, in comparison with Nambu-Goto type approaches. 

We begin our related considerations by expanding the action (71) around the helicoid 
classical solution: 

S = Sd--j dr y{T)x''\T) + ^ j dr j dfyar) x 



X 











1, d-'^ 5^A[C] 



(106) 



where y = x — x^^. 

Using the results of Sections 2 and 3 one can easily determine that 

S^A[C] 



SXa{T)Sxp{f) 



x,(r)i:.(f)A2,.^[x(r)-a;(f)]- 



§^S{r-f) J dS,,{z')A^:i,^[z{C'-x{T)] + 



S(C) 



+Xa{T) j dsa{f,s)zx{f,s)z'p{f,s)t''''^^'A^aMf,s)-x{T)], (107) 



where we have written 



Sx„ 



= Sf^Jir - f)a(f,s). 



The second term on the rhs of Eq.(107) is simply the area derivative which, as we have seen 
in Section 3, has the general form ~ gaX^ — ggXa- Thus, for the classical solution g[x'^^] 
it gives zero contribution. It is, furthermore, easy to verify that the third term in (107) also 
disappears for x — x'^^. We, therefore, conclude that 

5''A[C] 



X 



-^''T)x1^{f)A^Xp\.^'^\T)-x{f) 



(108) 



5xi{T)5x${f) 

Inserting Eq.(107) into Eq.(108) and taking into account that the dominant contribution 
to the two-point correlator comes from the region r ~ f we find 

b 



S ^ S^x + j daya{cr) 



1 mo 



2 I X 



A2 

^9 



(109) 



Let it be remarked that to arrive at the above relation we have adopted the expansion of the 
two-point correlator indicated in Eq. (28). We have also used the helicoid parametrization 



28 



observing, at the same time, that the eikonal hnes give null contribution. One further realizes 
that the contributions of the two helical curves to the linear term in (106) cancel each other, 
since x^^^^s) = -x'^^\s) and r2 - ri ~ r4 - ra ~ 

The non-trivial contribution of the helical curves is incorporated in the term 

^a, = S^,-^,{x^^xf+X^^X^;^), (110) 

the origin of which is the second functional derivative, cf. (106). The mass parameter in 
(109) has the same source and is defined as 



=1 i: I T^J dw (d{w') + D,{w') + ^D,{w')\ . (Ill) 
^ ^ 

The differential operator entering Eq.(109) has no zero eigenvalues since the "classical" 
solution is, in fact, the one that annihilates the g^-function. Accordingly, the calculation 
of the path integral over y = x — x'^^ does not require any particular regularization. A 
straightforward calculation shows that 

deta;«, = -^l--j = ^^^. (112) 

Thus the matrix LOa.p can be diagonalized and the y-integral can be easily performed. How- 
ever, in the limit mo ^ it can be immediately seen that the integration over the boundary 
fluctuations gives prefactors which are powers of the logarithm of the incoming energy and 
as far as Regge behavior is concerned, they cannot change the behavior that was determined 
in the previous section. 

The next task is to take up the issue of the spin-field dynamics contribution to the 
scattering amphtude. As seen in Section 3 a spin factor is associated with each segment of 
the worldline path. This factor receives contributions from two sources. The first one is 

/dr / dS. A('\z -x)-J^j drM^i^^^ l7„7.] (113) 
S{C) ^ 

and is obviously zero for the classical trajectory (72). 

The other term has the form 

P=\jdrj dr'J^^A^l^^ix - x')Jxp = U dr J dr\D + D^) + ^ j dr j dr'^x - x'fD[. 

(114) 
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In the stochastic hmit, within which we are working, the integrals in the above equation 



give appreciable contribution only for 



XiT 



X(T 



X 



t' |< To. More con- 



cretely, consider the contribution to (114) from the helical curve (1 — > 2). A straightforward 
calculation shows that the analytically continued result is 



p = -{t2-hy 



X 



(115) 



where we have written t — it for the time variable and denoted 



J dwD(w^) 

I 

'^5 

/ dww'^Di{w'^) 
\ 



1/2 

n \ 



y) 



(116) 



As has been mentioned in Section 3 and discussed in [27], contribution (115) has an 
interesting role as far as the form of the fermionic propagator is concerned. As it has been 
shown in [26,32] this non-perturbative "paramagnetic" contribution corrects the self-energy 
of a bound light quark and consequently the Regge intercept without changing the slope. 
We shall not discuss here this interesting issue leaving it for a forthcoming study. In this 
paper we bypass the problem using the parametrization [17,36]. 

The remaining spin structure is summarized in the chain 



i=4 



1 - 



2m, 



;ii7) 



which must be sandwiched between the external spinor wavefunctions representing the in- 
coming and outgoing quarks (in the simple picture wherein the meson wavefunction is just 

the product of free spinors). The non-trivial dynamics of the process are now incorporated 
into the fact that the vectors x^^\ i = 1,2,3,4 forming the boundary of the helicoids, are 
3-dimensional vectors with | x^'^^ \'^— const. For i — 1,3 turns the factor in (117) to the 



operator 1 — 



TP 



(i) 



\p(i)\ ■ 

For i — 2,4: the matrices 



1 - 



1 



and 



74 = 1 



2m T2 — Ti 

b 



x^^\b) 



1 - 



7-x(2)(6) 



X 



1 



2m T4 — T3 



j-x^^\b) ^ 1 



7-i;(^)(6) 



X 



(118) 



(119) 
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are also representations of projection operators. As shown in [31] the matrices (118) and 
(119) are the direct product of two 2x2 matrices each of which are by themselves projection 
operators. Given these observations it becomes a matter of simple algebra to find that the 
standard kinematics are reproduced. 

6. Concluding Remarks 

The central objective of this paper was to assess the merits of the Stochastic Vacuum 
Model of Dosch and Simonov as a credible representative of QCD. From a methodological 
standpoint we employed the path-integral approach for the casting of the theory, a practice 
that has been proved an ideal tool for the exploration of its non-perturbative aspects on 
which the present study is focused. In the first part of the paper-and at a purely theoretical 
level-we verified both the loop equations and the Bianchi identity through the SVM, an 
occurrence which further solidifies the credibility of the model. Wc have also derived an ex- 
plicit expression for the spin factor that represents the non-perturbative spin-field dynamics 
and necessarily enters the analysis of physical processes. In the second part of the paper 
we assessed the effectiveness of the SVM, always in its path-integral casting, to confront 
a dynamical problem where a non-perturbative treatment is essentially important. More 
explicitly, we appropriately modelled a meson-meson scattering process in the Regge kine- 
matical regime. In a "semiclassical" approximation and always working in the framework 
of the Stochastic Vacuum Model we found a Regge-type behavior for the scattering ampli- 
tude with linear Regge trajectories. The specific methodology we followed is entirely based 
on the capability of the SVM to represent the non-perturbative content of QCD and, per- 
haps, it traces a way for analytically calculating Regge trajectories in the physical region of 
scattering, i.e., square momentum transfer t < 0. 
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Appendix A 



We give here the proof of relation (20) which appears in the text and whose role is signifi- 
cant for the derivation of the equation that determines the surface on which the two-point 
connector "lives". We begin by writing the expression for the connector in Eq. (8): 

1 



(f){z, xq) — Pexp ( i j dw -.41= Pexp i j dTw{T) ■ A{w{t)) 



(A.l) 



where tyju(O) = Xq^, 'U'^(l) = z^. Taking the functional derivative of (A.l) we find 



or 



T^^^j dTPexp(^i J dTw-A{w)^[dr5{T—T')+5{T—T')w^dij_Ai,{w)]Pexp(^i J drw ■ A{w) 

' (A.2) 



- = iS{l-T')A^{z)Pexp 



z 

i j dw ■ A{w) 



L xo 



i5(r')Pexp 



i j dw ■ A{w) 

. xo 

j dw ■ A{w 



xo 



Sw^{t') 

+ igw,^{T')P exp [^i J dw ■ A{w)j F^,^{w{t'))P exp 

It, accordingly, follows that the variation of the connector reads 

Sep = A^{z)5z^,(t){z, ,To) - i(f){z, xo)Af,ixo)5xo^^ + 

+ [ dTW,,{T)4){z,w{T))F^^{w{T))4){w{T),XQ)5w^{T). 

Jo 

Keeping everything but the end point constant one immediately deduces that 



Af,{xo) + 
(A.3) 



(A.4) 



dz„ 



dw\ 

Ai,{z)4>{z,xq) - dTW^{T)4>{z,w{T))F^x{w{T))4>{w{T),XQ)- — . (A. 5) 

Jo oz,. 
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Appendix B 



In this Appendix we present a parametrization of the functions D and Di used extensively 
in the present paper. This parametrization is supported by lattice data and is extensively 
discussed in Refs. [17,36]. 

The exact relations defining the functions are 

7:\Nl-l) G, 7:\Nl-l) G, 
^ = 2Nc 2Nc ^^l-'^)^^'^' ^^-1) 

where Djv and Di^at are functions which determine the structure of the two-point correlators, 
as defined in [17]. The factor G2 is defined as follows 

G, ^ (0 I £-,F^MF^M I 0) = ^A(f),,(0). (B.2) 

For the above correlator we shall adopt the value given in Ref [8] , namely G2 — {0A96)'^GeV^. 
The value of the numerical quantity k, in (B.l) is estimated in the same reference to be 0.74. 
The ansatz for the function is [17] 

27 1 



d4 a"^ J 



4 ' 



(B.3) 



where 



00 



a = J dzDN{z). (B.4) 


A simple calculation shows that 

1 3tt 
Dn{z) ^ wKi{w) - -w'^Ko{w), w = — \z\, (B.5) 
4 8a 

with Ki, denoting a Bessel function. The correlation length Tg can be deduced from Eq. 
(B.5): 

^ 8a s 

r, = -. (B.6) 

The estimated value of a is 

a f« 0.35fmor Tg fti 0.297fm. (B.7) 
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With the help of ansatz (B.3) and using (B.7) one can determine the string tension: 



or 



^ 2Nc 24 

1^2^'(2iV^-l)G2 



wKi{w) — -w^Ko{w) 
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-— -K27r = a^GsKTT— 0.175GeV^ 

2Nc 24 81 



The anzatz for the function i^i,ivis deduced from the equation [17,36] 

d 



4 + ^. 



or 



Di,n{z) = \ r dw[Aw^K,{w)-w'Koiw)] 
z Jo 



The coefficient of the rigidity term entering Eq. (94) can now be calculated: 

- = ^T^d'ww'[2D,{w)-D{w)]^ 
an 62 J 

J d^ww'^[2{l - k)Di^n{w) - kDn{w)] 



32 3 2Nn 24 



= Irf '^g-"gi2(l - .)32. « 0.276. 

32 ^ 2Arc 24 ^ ' 
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